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One of the most promising probes to study deconfined matter created in high energy nuclear 
collisions is the energy loss of (heavy) quarks. It has been shown in experiments at the Relativistic 
Heavy Ion Collider (RHIC) that even charm and bottom quarks, despite their high mass, experience 
a remarkable medium suppression in the Quark Gluon Plasma. 

In this exploratory investigation we study the energy loss of heavy quarks in high multiplicity proton- 
proton collisions at LHC energies. Although the colliding systems are smaller than compared to 
those at RHIC (p+p vs. Au+Au) the higher energy might lead to multiplicities comparable to 
Cu+Cu collisions at RHIC. Recently the CMS collaboration has shown that these particles likely 
interact among each other. The interaction of charm quarks with this environment gives rise to a 
non- negligible suppression of high momentum heavy quarks in elementary collisions. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Nq 
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The start of the Large Hadron Collider (LHC) gives ac- 
cess to an unexplored area of high energy physics. While 
it has been commonly assumed that only high energy 
nucleus- nucleus collisions are able to produce a decon- 
fined state of nuclear matter, commonly acknowledged 
as the Quark Gluon Plasma (QGP), the high energy and 
the resulting high multiplicities in proton-proton colli- 
sions at the LHC might allow for a QGP as well. The 
multiplicity of 7-14 TeV proton-proton collisions are com- 
parable to Cu-|-Cu collisions at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL). Since one expects at least a partially deconfined 
state in Cu-|-Cu collisions at an center-of-mass energy of 
200 AGeV [ll-@] it might be possible to reach a phase of 
deconfined matter in p-l-p collisions at the LHC as well. 
Recently the CMS collaboration indicated this possibil- 
ity experimentally Q. However, the proton-proton data 
always served as a baseline during the RHIC program 
and thus it might prove hard to find proper observables. 

Heavy mesons would be an ideal probe to verify the 
existence of a QGP in proton proton collisions. If a de- 
confined plasma is created and the heavy quarks pass 
this plasma, the same strong interaction between plasma 
constituents and heavy quarks is expected which one has 
observed in RHIC heavy ion data. If no plasma is cre- 
ated and all hadrons are created as color-neutral parti- 
cles in the vacuum no modification is expected. This 
argument originates from the fact that the interaction of 
the heavy mesons with the surrounding hadrons is much 
weaker compared to the strong interaction heavy quarks 
undergo in a deconfined state of matter. Thus a clear 
distinction for experimental signals is given. 

In this exploratory study we show that in very high 
multiplicity proton-proton collisions one can measure 
a small, however non- negligible energy loss of heavy 
quarks. For our studies we use the model proposed in 
|8l4l(J| and adjust it to proton-proton collisions at LHC. 
The purpose of this letter is to demonstrate the effect 
of heavy quark energy loss in proton-proton collisions at 



LHC energies. It is organized as follows: First, we will 
describe the models we use to calculate the underlying 
medium properties from the proton-proton collisions. We 
will then present the selection of high multiplicity proton- 
proton event within our theoretical framework. After 
that, we briefly review the elementary reactions that are 
the basis of our energy loss calculation. Finally we dis- 
cuss the energy loss of heavy quarks in proton-proton 
collisions and discuss possible experimental observables. 
The letter ends with conclusions. 

For our analysis, we split the medium description from 
the description of the elementary interaction of the heavy 
quarks with the medium (i.e. the light quarks and the 
gluons). Finally we combine both and calculate the spec- 
tra and a quantity similar to the nuclear suppression fac- 
tor Raa in heavy ion collisions to specify the medium 
modification. 

In order to simulate the proton-proton collisions and the 
medium that is created in those collisions we utilize the 
EPOS model including a hydro-dynamical phase, intro- 
duced in e.g. flll - fisj . It has been thoroughly tested for 
proton-proton collision physics and is known to describe 
p+p collision data within reasonable bounds. 

EPOS is a consistent quantum mechanical multiple 
scattering approach based on partons and strings |14| . 
where cross sections and the particle production are cal- 
culated consistently, taking into account energy conserva- 
tion in both cases (unlike other models where energy con- 
servation is not considered for cross section calculations 
[iBl)- A special feature is the explicit treatment of projec- 
tile and target remnants, leading to a very good descrip- 
tion of baryon and antibaryon production as measured in 
proton-proton collisions at 158 GeV at CERN [ll|. Nu- 
clear effects related to CRONIN transverse momentum 
broadening, parton saturation, and screening have been 
introduced into EPOS [l7[. In heavy ion collisions (and 



more recently also in proton-proton collisions) collective 
behavior is taken into account [l^ , in the following fash- 
ion: the initial scatterings, as described above, lead to the 
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formation of strings, which break into segments, usually 
identified with hadrons. When it comes to heavy ion col- 
lisions, the procedure is modified: one considers the situ- 
ation at an early proper time tq, long before the hadrons 
are formed: one distinguishes between string segments in 
dense areas (more than some critical density po segments 
per unit volume), from those in low density areas. The 
high density areas are referred to as core, the low density 
areas as corona [l^. It is important to note that initial 
conditions from EPOS arc based on strings, providing 
a flux-tube like structure in case of individual events (a 
single flux tube in case of many overlaid events, for a 
schematic view of a fluxtube we refer to Fig. [Ij. Based 
on the four-momenta of the string segments which consti- 
tute the core, we compute the energy density £{tq,x) and 
the flow velocity v{to,x). Having fixed the initial condi- 
tions, the system evolves according to the equations of 
ideal hydrodynamics. 
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FIG. 1. (color online) Schematic view of a parton ladder. 

In order to identify high multiplicity proton-proton col- 
lisions we trigger on the number of elementary collisions 
(labeled v) within the proton-proton collision itself, i.e. 
the number of strings or color fluxtubes formed in the col- 
lisions. A fixed number of elementary collisions v implies 
a distribution in multiplicity (as well as in impact pa- 
rameter) which has a certain width. The distribution of 
impact parameters for two different choices of the number 
of binary collisions is shown in Fig. [2] While this is not 
directly accessible within an experimental setup, an aver- 
age multiplicity can be calculated and thus compared to 
experiment. For v = 5 the average multiplicity of charged 
particles at mid-rapidity dNch/<iy\(y^Q-j is 16.05, for col- 
lisions with 1/ = 10 it equals 29. We are calculating the 
energy loss for matter with an energy density larger than 
1.5 GeV/fm^ Between 1.5 GeV/fm^ and 0.4 GcV/fm^ 
we are treating the matter as a so called mixed phase, 
with a reduced rate for reactions of heavy quarks with 
the medium. These values are in accordance with former 
hydrodynamical calculations applied to Au-f An reaction 
at RHIC. The impact of the hadronic phase has been 
checked and can be safely neglected. The interesting, 
possibly deconfined, phase is then roughly 2-3 fm in di- 



ameter (see e.g. ^^■) It has been cross-checked that for 
smaller values of the number of elementary interactions 
(i.e. smaller on the order of 1-2, which corresponds 
to multiplicities below 10) the medium modification of 
heavy quarks vanishes. A discussion on experimental ob- 
servables can be found at the end of this letter. 
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FIG. 2. (color online) Impact parameter distribution for two 
different choices of tiie number of elementary scatterings (flux- 
tubes). 

The energy loss model we use for our studies is 
MC@sHQ, for details we refer to For brevity we 

do only repeat the key features of the model. Except 
for adjusting aspects of the technical implementation the 
physical picture remains the same in p-t-p collisions as in 
Au-|-Au collisions. 

The two major improvements in comparison to conven- 
tional pQCD inspired models are the inclusion of a run- 
ning coupling constant and a proper matching of the in- 
frared regulator of the gluon propagator to HTL calcula- 
tions. These changes allowed to describe data in heavy 
ion collisions at RHIC within reasonable K-factors. For a 
recent comparison to experimental data, we refer to (lo| . 
Additionally to coUisional energy loss a radiative compo- 
nent to the energy loss is implemented. This is calculated 
as incoherent gluon bremsstrahlung pushing to its limit 
the idea that the gluon formation time is strongly re- 
duced when the gluon is radiated off heavy quarks. For 
more information we refer to [13 . The probability of col- 
lisional or radiative energy lost per unit length for a 10 
GcV c-quark is shown in Fig. [3] as a function of energy. 
CoUisional energy loss (dashed line) mostly consists of a 
large number of collisions in which the incoming heavy 
quark loses a small part of its energy, hence the pro- 
nounced peak at around zero in Fig. [S] For radiative 
losses (solid line) the probability to lose a larger fraction 
of the heavy quark energy in a single process is more 
important. 

To combine the description of the medium and the 
energy loss, we generate an averaged medium, i.e. an 
energy density profile with EPOS. Within this medium 
heavy quarks are generated and then propagated accord- 
ing to the Boltzmann equation. We then compare two 
different scenarios. First, the heavy quarks lose energy 
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according to the description in [T^]- Alternatively the 
heavy quarks move freely in the medium. The ratio of 
both final spectra - that is, the quenching - will then 
provide a measure for the lost energy. 
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FIG. 3. (color online) Probability of coUisional (dashed line) 
and radiative (solid line) energy loss per unit length as a func- 
tion of energy loss uj. 
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FIG. 4. (color online) The Raa equivalent in proton-proton 
collisions for different choices of the number of elementary 
scatterings (fluxtubes). For v = h the average multiplicity 
of charged particles at mid-rapidity dNch/dy\{y=o) is f6.05, 
for collisions with v = 10 it equals 29. The curve la- 
beled Rhm/lm depicts a high multiplicity spectrum = 10, 
dAfch/dj/|(j,=o)=29) divided by a low multiplicity spectrum 
{^=1, dNch/<iy\(y=o) ~ 3.62). Error bars shown are statisti- 
cal. 



In Fig. |4] we show the quenching (the equivalent of 
a nuclear suppression factor Raa for proton-proton col- 
lisions) for charm quarks. Theoretically this is well- 
defined, since one can calculate the medium-influenced 
spectra as well as the vacuum spectra by neglecting rein- 
teraction of the heavy quark with the medium. The ra- 
tio of both then gives the desired quantity. Additionally 
the below defined observable Rhm/lm is shown, which 
shows a very similar suppression. One should note that 
the suppression for v = b and = 10 seem very simi- 
lar. This originates from the fact that the heavy quark 
undergoes few interactions due to the generally small in- 
teraction zone. Once the energy density reaches a certain 
threshold the suppression sets in and the particles with 
high transverse momentum escape the high energy den- 
sity zone rapidly. Thus the low pT part of the spectrum 
is effected more, which is in line with the observation in 

Fig. m 

One observes a small, however non-negligible energy 
loss of charm quarks at high transverse momentum 
manifesting itself as a 10% suppression of the spectra at 
large pt {pt > 7GeV). This prediction should however 
be considered with a grain of salt: it is known - for a 
recent review, see [l^ - that the heavy quarks produced 
inside a hot medium undergo a transient regime before 
losing energy proportional to the path-length, an effect 
not implemented in our simulations. Also additional 
initial state effects might have shift the data in the 
opposite direction and more systematic studies are 
possible once LHC data arrives. 

When trying to define this observable for experimen- 
tal applications, we are confronted with the complication 
that the reference measurement uses the same colliding 
systems as the medium measurement. Thus we need a 
quantity which refers to experimental observables. We 



suggest a new ratio, which can be measured experimen- 
tally. It requires knowledge of multiplicity of the proton- 
proton collision, and a clear cut between low multiplicity 
and high multiplicity collisions. Then one can construct 
an observable similar to Rc p in heavy ion collisions (20| . 
Measuring the spectra in high multiplicity p-t-p collisions 
and dividing them by the spectra obtained from low mul- 
tiplicity collisions one constructs 

dN/dpT,HM 

r? Nhm {-\ \ 

Rhm/lm - dN/dpT,LM " 

Nlm 

Here, dN/dpr.HM denotes the HQ spectrum ob- 
tained from high multiplicity proton-proton collisions, 
dN/dpT^LM is the spectrum as obtained from low multi- 
plicity proton-proton collisions. Nhm and Nlm denote 
the multiplicities in high multiplicity or low multiplicity 
proton-proton collisions respectively. In our analysis 
we assume the initial HQ spectrum scales linearly with 
the number of binary reactions. The advantage of this 
observable is that it is " self- normalizing" , due to the 
scaling with the respective multiplicities, which serve 
as a normalization of the spectra. In principle one 
should aim for the bins with the highest and lowest 
multiplicity, for our analysis we have chosen two values 
diVch/dy|(j^=o) = 29 {i^ = 10) and 16.05 = 5) in 
the high multiplicity range to compare possible effects. 
For the low multiplicity bin we selected the one with 
'iNch/<iy\{y=o) = 3.62 {v = 1), however experimental 
statistics might be enough to aim for even higher 
multiplicities. Apart from experimental difficulties we 
predict a quenching of roughly 10% at pT > 10 GeV in 
high multiplicity proton-proton collisions as indicated in 

Fig. m 
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which can be evaluated in different multiphcity classes 
and thus different values of z. The scaling with z might 
give some insight into the centrality dependence of the 
energy loss of heavy quarks. 

Shown in Fig. [S]is a calculation of R(z) within our ap- 
proach. One observes a decrease of roughly 10-15% with 
increasing centrality/multiplicity of the collision, which 
is in line with the expectation from the previous argu- 
ments. 



CONCLUSIONS 



FIG. 5. (color online) High transverse momentum suppression 
of charm quarks as a function of 2: = Nch/ {Nch}- Shown is 
the medium calculation divided by the vacuum calculation, 
averaged above 10 GeV. Error bars shown are statistical. 



Instead of a low multiplicity sample one could choose 
minimum bias data as well, since it is dominated by col- 
lisions producing few number of particles. Although this 
might be more accessible to experiments, we suggest that 
a clear cut between low and high multiplicity collisions 
enhances the viability of this observable. 

Another possible scaling is a scaling with z = 
Nch/ (Nch), where Nch denotes the number of charged 
particles in a certain centrality class and (Nch) the mean 
value from minimum bias (MB) data. A possible experi- 
mental observable would be 



In our exploratory study we have shown that if in high 
multiplicity proton-proton collisions a deconfined phase 
is produced, heavy quarks will be a suitable probe for its 
experimental analysis. If fast heavy quarks pass a de- 
confined medium they lose energy due to the interaction 
with the medium and thus are quenched. This quench- 
ing reaches about 10% at large transverse momentum. 
By dividing the heavy quark spectra for high and low 
multiplicity p-l-p collisions (scaled accordingly) this ef- 
fect can be experimentally observed at LHC experiments. 
An insightful measurement would be the energy loss as a 
function of multiplicity, e.g. as a function of the variable 
z ~ Nch/ {Nch)- By this measurement a clear onset of 
medium formation in proton-proton collisions could be 
observed. 
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